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ABSTRACT 

In close eclipsing binaries, measurements of the variations in binary's eclipse tim- 
ing may be used to infer information about the existence of circumbinary objects. To 
determine the possibility of the detection of such variations with CoRoT and Kepler 
space telescopes, we have carried out an extensive study of the dynamics of a binary 
star system with a circumbinary planet, and calculated its eclipse timing variations 
(ETV) for different values of the mass-ratio and orbital elements of the binary and the 
perturbing body. Here, we present the results of our study and assess the detectability 
of the planet by comparing the resulting values of ETVs with the temporal sensitiv- 
ity of CoRoT and Kepler. Results point to extended regions in the parameter-space 
where the perturbation of a planet may become large enough to create measurable 
variations in the eclipse timing of the secondary star. Many of these variations point 
to potentially detectable ETVs and the possible existence of Jovian-type planets. 

Key words: methods: numerical, techniques: photometric, binaries: eclipsing, plan- 
etary systems, planets: detection 



1 INTRODUCTION 

Approximately 70 percent of the main and pre-main se- 
quence stars are members of binary or multiple star sys- 
tems. Observational evidence indicates that many of these 
systems contain potentially planet-forming circumstellar or 
circumbinary disks implying that planet formation may 
be a common phenomenon in and around binary stars 
jMathieu lll994 !:Ukes on et al. lll998l:lRodriguez et al. 1 11998: 
White et al"ll 999: Sil bert et al. ll2000l : lMathieu et al. IboOC ; 



TriUing et ani2007l ) 



Many efforts have been made to detect such binary- 
planetary systems. In the past two decade, even though sin- 
gle stars were routinely prioritized in search for extrasolar 
planets, many of these efforts were successful and resulted 
in the detection of approximately 40 planet-hosting binary 
star systems. The discovery of these systems have led to 
the speculation that many more planets may exist in and 
around binaries and prompted astronomers to explore the 
possibility of the detection of these planets with different 
detection methods. We refer the reader to Planets in Bi- 
nary Star Systems (|Haghighipour I [20101 ^ for an up to date 
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and comprehensive review of the current state of observa- 
tional and theoretical research in this area. 

In general, one can distinguish three types of planetary 
orbits in a binary star system (Fig. [l]): 

(i) S-Type, where the planet orbits one of the two stars, 

(ii) P-Type, where the planet orbits the entire binary, 

(iii) T-Type, where the planet orbits close to one of the 
two equilibrium points L4 and L5 (Trojan planets). 

Currently, the most known planets in binary sys- 
tems are in S-type or bit s [for mo re details see e.j 
'Pilat -Lohinger fc DvorakI (|2002| '1. iPilat-Lohinger et al.l 



(I2OO3I ), and iHaghighipour et al.l ((20101)]. The stellar sepa- 
rations of many of these binaries are larger than 100 AU 
implying that the perturbation of their farther companions 
on the formation and dynamical evolution of planets around 
their planet-hosting stars may be negligible. However, in 
the past few years, radial velocity and Astrometry surveys 
have been able to identify five binary star systems with 
separations of approximately 20 AU where one of the stars 
is host to a Jovian-type planet. These syste ms, namely, 
GL 86 (iQueloz et al. I I2OO0I : IEIs et al. I I2001I). y Cephe i 



( Hatzes et al. I |2003|) 



iRaghavan et al. I l2006fl 



HD 
HD 



41004 
196885 



(iZucker et al. 



and HD 176051 l|Muterspaugh et al.l [20T0F present unique 



Correia et al. 



20041 : 



liooi). 
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Figure 1. Schematic illustration of S-type and P-type orbits in a 
binary star system. The yellow and red circles represent the stars 
of the binary revolving around their common center of mass (the 
blue circle). 



cases for the study of planet formation and dynamics in bi- 
naries as the perturbation of their secondary stars will have 
significant effects on the dynamics of their circumprimary 
disks and their capability in forming planets. 

The success of radial velocity and astrometry in detect- 
ing planets in S-type orbits raises the question that whether 
other techniques can also detect planets in and around bi- 
nary stars. A recent success has come in the form of the 
detection of the first P-type planets using the ETV method. 
Modelin g the variations in the e clipse timing of the binary 
NN Ser. iBeuermann et al.l l|2010l ) have shown that this sys- 
tem is host to two planets with minimum masses of 2.2 and 
6.9 Jupiter-masses in a 2:1 mean-motion resonance in cir- 
cumbinary orbits. 

Our work has been motivated by this discovery. NN 
Ser (ab) is an eclipsing short period binary which shows 
long-term ETVs. The detection of P-type planets around 
this system suggests that close, eclipsing binaries may be 
promising candidates in the search for new exoplanets. As 
many of these bin aries lie in t he discovery space of CoRoT 
l|Goldmannll2010h and Kepler l|Coughlin et al.ll201(]| ) space 
telescopes, we focus our study on exploring the prospects of 
the detection of such objects. 

The idea of photometric detection of extrasolar plan- 
ets around eclipsing binarie s was first presented by 
ISchneider fc ChevretonI (|l990l ). This idea that was later 
developed by rn a ny authors 



Dovle et al.l Il998l: 



Muterspaugh et al.l 



1995 



2008 



Dovle &: Dees 

is based on the fact that a cir- 



cumbinary planet can perturb the orbit of the two stars and 
create variations in their eclipse timing. The measurements 
of these variations, when compared with theoretical models, 
can reveal information about the mass and orbital elements 

of the pertur bing planet. 

Recently ISvbilskv et al.l (|2010l ) studied the potential of 
the ETV method in detecting giant planets in circular, cir- 
cumbinary orbits. These authors have shown that ground- 
based photometry may have advantage over CoRoT and 
Kepler space telescopes in detecting circumbinary planets. 
They suggest that the selective nature of the target stars in 
the fields of view of CoRoT and Kepler causes the ETV sig- 
nals of eclipsing binaries in these regions to be smaller than 



the detection sensitivities of these telescopes. These authors 
also suggest that if target binary stars include presumably 
less stable contact binaries, the CoRoT's and Kepler's capa- 
bilities of detecting of circumbinary planets increase by four 

times. 

In this paper, we extend the study by ISvbilskv et al.l 

l|2010l ) to binaries with planets in eccentric and reso- 
nant orbits. It is expected that similar to the transit 
timing variation method where TTV signals are strongly 
enhanced when t he transiting and perturbing planets 
are in resonance jAgol et all l2005l : lAgol fc SteffenI l2007l : 
ISteffen fc Agoll |2005| ), a resonant perturbing circumbinary 
planet will also produce high ETV signals. The goal of our 
study is to identify regions of the parameter-space for which 
this signal is within the temporal sensitivity of CoRoT and 
Kepler space telescopes. 

The outline of our paper is as follows. We present our 
models in section 2 and discuss their stability in section 3. 
Section 4 has to do with the calculations of ETVs and the 
prospects of the detection of P-type planets. Section 5 con- 
cludes this study by summarizing the results and discussing 
their implications. 



2 MODEL 

We consider a close, eclipsing binary with a giant planet in 
a P-type orbit. To ensure that the effect of the variations 
of the mass-ratio of the binary is included in our study, we 
consider the following three models where mi and are 
the masses of the primary and secondary stars, respectively, 

• model 1: mi = m2 — 0.3 Msun, 

• model 2: mi = l,m2 — 0.5 Msun, 

• model 3: mi = m2 = 1 Msun- 

As most of the stars in the solar neighborhood are of spectral 
type M, this choice of models ensures that low-mass binary 
stars are also included in our simulations. 

Eclipsing binaries are morphologically classified as 

(i) detached systems, if neither component fills its 
Roch lobe (separated stars), 

(ii) semi-detached systems, if only one component fills 
its Roche lobe, and 

(iii) over contact systems, if both components exceed 
their Roche lobes. 

The Roche lobe marks the volume limit at which the star 
may begin to lose substantial amount of matter to its com- 
panion. In this study we consider a detached binary with 
an initial separation of 0.05 AU. For the binary models 1 and 
3, this separation corresponds to a perio d of approx i matel y 
5 and 3 days, respectively. As shown by iGoldmannI (|2010t ). 
most candidate eclipsing binaries in the discovery space of 
CoRoT have periods between 1 and 10 days. 

When dealing with close binaries such as the models 
considered here, the intrinsic eccentricity of the binary (cbin) 
can be neglected. This is due to the high probability of cir- 
cularization that is caused by interstellar tidal forces. In this 
study, we consider the initial orbital eccentricity of the bi- 
nary to be zero or have very low values. The latter is due 
to the fact that the timescale of c ircularization is depend ent 
upon specific stellar parameters (|Zahn fc Bouchetll 19891 ). 
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Figure 2. The emax maps of a circumbinary planet in the binary model 1 (loft column), 2 (middle column), and 3 (right column). The 
top row shows the maps for a 1 Mj planet, the middle row is for a 5 Mj planet, and the bottom row corresponds to a 10 Mj object. 
The violet and blue regions represent small values of emax and correspond to bound orbits whereas the yellow region points to escape. 



3 STABILITY ANALYSIS 

Prior to calculating the variations in the eclipse timing, it is 
necessary to identify regions where a P-type planet can have 
a stable orbit. For this purpose, we integrated the Newtonian 
three-body system of the binary and its planet for different 
values of the planet's mass and orbital elements. The mass 
of the planet was taken to be 1, 5, and 10 Jupiter-masses. 
The values of its semimajor axis and orbital eccentricity with 
respect to the barycenter of the system were varied between 
a — 0.1 — 0.3 AU and e = — 0.3 in increments of Aa — 0.01 
AU and Ae = 0.01, respectively. 

Numerical integrations were carried out u sing the Lie- 
series metho d and Bul irsch-Stoer integrator (Lichtoncggei' 
[l984; Hansl meier fc Dvor ak 1985.; Eggl fc Dvor ak 2010). We 
integrated the system for lO'' years which corresponds to the 
following number of periods of the planet at 0.1 AU: 

• model 1: 2.45 x 10^ periods, 

• model 2: 3.90 x 10^ periods, 

• model 3: 4.50 x 10^ periods. 

The stability of the planetary motion was controlled 
by measuring the value of its maximum orbital eccentricity 
{emax)- We monitored the changes in the eccentricity of the 
planet throughout the integration and determined its high- 



est value. If the orbit of the planet became parabolic (e = 1), 
we considered the system to be unstable. 

Figure 2 shows the values of planet's emax at different 
distances from the barycenter of the binary. The left column 
corresponds to the binary model 1, the middle column to 
binary model 2, and the right column to binary model 3. 
Also, from top to bottom, each row corresponds to the values 
of emax for a 1 Mj, 5 Mj, and 10 Mj, respectively. In each 
panel, the axes represent the initial values of the semimajor 
axis and eccentricity of the planet at the beginning of the 
integration. The color at each point depicts the maximum 
value that the eccentricity of the planet acquired during the 
integration. As indicated by the scale on the right side of the 
figure, yellow corresponds to parabolic orbits and denotes 
instability. 

An inspection of the emax maps shown in Fig. 2 indi- 
cates that for planets in circular or low eccentricity orbits, 
the range of the semimajor axis for which the planetary or- 
bit may be bound is large and does not change much for 
different values of the planet's mass. For instance, as shown 
by the top left panel of this figure, the boundary of the 
unstable zone for a 1 Mj planet in a circular orbit in the 
binary model 1 is interior to 0.1 AU. This suggests that all 
circular orbits with semimajor axes larger than 0.1 AU in 
this model may be bound. As the mass of the planet grows 
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(middle and bottom panels of the left column in Fig. 2), 
the outer boundary of the unstable region slowly progresses 
toward larger semimajor axes. Such a trend is also seen in 
the emax maps of binary models 2 and 3. As shown by the 
lower panels in the middle and right columns, the bound- 
ary of the unstable zone for circular orbits is slightly shifted 
outward to 0.11 AU. These results are consistent with the 
results of the stability analysis of a test particle in a cir- 
cumbinary orbit as presented by Dvorak et al. (1989) and 
Holman & Wiegert (1999). According to the estimate of the 
boundary of the stable zone as given by these authors, the 
critical distance beyond which the orbit of a Jovian planet 
in all our three binary models will be stable is in the range 
of 0.1 to 0.12 AU. 

While the value of emax can be used to identify parabolic 
(unstable) orbits, it cannot be used as a rigorous indica- 
tor of planet's orbital stability. A low value of Cmax im- 
plies a bound planetary orbit. But that is only for the du- 
ration of the integration, and there is no guarantee that 
the orbit of the planet will stay bound for a long time. In 
order to determine the orbital stability of the pl anet, we 
used the Fast Lyapunov Indicator (FLI) (Froesc hle et al.l 
1 19971 ). The FLI is a chaos indicator that measures the 
exponential divergence of nearby trajectories and distin- 
guishes between regular and chaotic moti on. For details 



of thi s technique, we refer th e re ader to [Fro eschle et al 
(| 19971 ) . iLeea fc Froeschli l|200ir i and lFouchardet al.. (2002) 
Applications of FLI to planetary motion in binary sta r 



2002 ). 



20101 1 



systems can be found in i Pilat-L ohinger fc DvorakI ( 
IPilat-Lohinger et al.l (l2003h andlllaghighi pour et al.l I 

Figure [3] shows a sample of an FLI stability map for a 1 
Mj planet in binary model 1. The colors in this figure rep- 
resent the values of FLI (the logarithmic scale on the right) 
which depict the degree of chaos. Dark colors correspond 
to less chaotic and regular orbits. As shown by this figure, 
there is a region in the parameter-space where the orbit of 
the planet is likely stable (dark area between 0.1 and 0.3 
AU). 

A comparison between the dark region of Fig. 3 and 
its corresponding region in the system's emax map in Fig. 2 
(upper left panel) points to an interesting observation: even 
though the maximum values of the planet's orbital eccentric- 
ity can reach to high values, the orbit of the planet may still 
be stable. In other words, elevated planetary eccentricities 
do not necessarily correspond to instability. For instance, as 
shown by the upper left panel of Fig. 2, orbits with initial 
eccentricities of 0.3 may have emax values up to 0.7. Never- 
theless, according to the FLI results shown in Fig. [31 they 
are probably long-term stable. 

Similar comparison can also be made between the 
chaotic region of Fig. 3 (indicated by red and yellow col- 
ors) and the system's emax map in Fig. 2. The results point 
to an inverse phenomenon, that is, instability among seem- 
ingly bound orbits. For instance, while for semimajor axes 
ranging from a = 0.11 AU to a = 0.14 AU, and the ec- 
centricities between e — 0.05 and e — 0.15, the Cmax map 
suggests bound planetary orbits, the FLI stability map indi- 
cates chaotic motion. The different color shades for chaotic 
orbits in the FLI map have resulted from the fact that the 
degree of divergence depends on the initial conditions of the 
orbits. Even though this shading may suggest less chaotic 
'islands' in parameter-space, when compared to the corre- 
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Figure 3. The FLI stability map for a 1 Mj P-type planet in the 
binary model 1. The violet and blue regions (small values of the 
FLI) correspond to regions of regular motion, and red and yellow 
denote chaotic orbits. 



sponding emax map, one can see that orbits inside the chaotic 
border zone become parabolic. 

Using the combination of emax and FLI analyses, we 
identified the stable planetary orbits in all our three binary 
models. As an example, the number of stable circular or- 
bits are shown in Table [1] As expected, stability is almost 
independent of the mass-ratio of the binary. 

We also carried out simulations for non-zero values of 
the binary's eccentricity. Figure 4 shows the maps of emax for 
a 1 Mj planet in the binary model 1. The binary eccentricity 
was chosen to be ebin ~ 0.05, 0.1, and 0.15. Table [2l shows 
the number of stable orbits for these simulations. As shown 
here, the unstable region expands out by more than two folds 
when ebin = 0.1. Our analysis indicated that for the value 
of the binary eccentricity ebin ~ 0.2, only 4 orbits remained 
stable (Table 2). 

The results of our stability analysis allow us to focus 
our calculations of ETVs on the region of the parameter- 
space where P-type planets have stable orbits. As shown, 
the planet's stable region shrinks for large values of its ec- 
centricity as well as the eccentricity of the binary. Initial 
binary eccentricities beyond ebin = 0.15 move this region 
to large distances (a = 0.24 AU) where the perturbing ef- 
fect of the planet on the dynamics of the binary becomes 
negligible (Fig. [3}. The closest possible stable orbit for a 
planet (a — 0.1 AU) is when the planet is Jupiter-mass, and 
both the planet and binary have circular orbits. As a result, 
for the purpose of calculating ETVs, we only consider fully 
circularized binaries. 



4 CALCULATION OF ECLIPSE TIMING 
VARIATIONS 

As indicated by the results of the stability analysis, a Jovian 
planet in a P-type orbit may be stable in the vicinity of a 
binary star system. Although small, the gravitational per- 
turbation of this planet may affect the motions of the two 
stars and cause their orbits to deviate from Keplerian. In an 
eclipsing binary, these deviations result in variations in the 
time and duration of the eclipse. Similar to the variations in 
the trans it timing of a planet due to a second perturber 
(Miralda -Escude fc Adams! [20051 : [Holman fc Murravl [20051: 
lAgol fc Steffen[ [2007[ : [Kippind [2009a[ [bl: [Ford fc HolmanI 
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Table 2. Number of stable orbits for binary model 1 and for 
different values of e^in- 
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Figure 4. The Cmax maps of a 1 Mj planet in the binary model 
1 with Bbin = 0.05 (top), 0.1 (middle), and 0.15 (bottom). 



Table 1. Number of stable orbits of all three binary models for 
Ebin = 0. 



Model 


A^planet 

{Mj) 


stable orbits 
(number) 


stable orbits 
(%) 


model 1 


1 


540 


82.9 




5 


530 


81.4 




10 


527 


80.9 


model 2 


1 


531 


81.5 




5 


529 


81.2 




10 


526 


80.7 


model 3 


1 


539 


82.6 




5 


538 


82.6 




10 


533 


81.8 



ebin 


A^planct 
(Mj) 


stable orbits 
(numbers) 


stable orbits 
(%) 


0.00 


1 


540 


82.9 


0.05 


1 


407 


62.5 


0.10 


1 


270 


41.4 


0.15 


1 


127 


19.5 


0.20 


1 


4 


0.6 



the variations in eclipse timing can be used to in- 
fer information about the mass and orbital elements of the 
circumbinary planet. In this section, we calculate the eclipse 
timing variations of our model binary star systems for differ- 
ent values of the mass and orbital parameters of the binary 
and its P-type planet. Our goal is to identify a range of 
these parameters for which the magnitude of ETVs will be 
within the temporal sensitivity of CoRoT and Kepler space 
telescopes. 

We simulated the dynamics of our binary models and 
their P-type planets in a barycentric coordinate system. 
Simulations were carried out for 1 year corresponding to ap- 
proximately 70 transits of the secondary star. We calculated 
ETVs by determining the difference between the eclipse tim- 
ing (ti) of the unperturbed system (star-star) and its corre- 
sponding value {t2) in the perturbed case (star-star-planet j3- 
Additionally, we subtracted the planetary induced constant 
rate of apsidal precession. We also Fourier-analyzed the re- 
sulting signal and compared the superposition of the main 
frequencies to simple estimates on maximum amplitudes 
[{dtmax — dtmin)/^] for the entire duration of integration. 
In cases where the differences between Fourier-composition- 
derived ETVs and maximum amplitudes were lower than 
10%, we chose to present the ETVs' maximum amplitudes. 



4.1 Implication for the detection of circumbinary 
planets 

As mentioned earlier, we would like to identify regions of the 
parameter-space for which the ETV signals of an eclipsing 
binary will be detectable by CoRoT and Kepler space tele- 
scopes. Since as indicated by the stability analysis, stable 
planetary orbits move to large distances as the eccentricity 
of the binary increases, we limited our calculations to circu- 
lar systems. We recall that in all our models, the separation 
of the binary is 0.05 AU. 

For the planet's motion, we considered both resonant 
and non-resonant orbits. Similar to transit timing variation, 
ETV signals are strongly enhanced when the P-type planet 
and the binary are in a mean-motion resonance (MMR). 
Results of our stability analysis indicate that the location of 
the 2:1 MMR at 0.079 AU is too close to the binary's center 
of mass to be stable (Dvorak et al. 1989, Holman & Wiegert 
1999). We, therefore, considered cases where the planet is in 
a 3:1 (0.104 AU) and a 4:1 (0.126 AU) MMR. 



^ We chose to integrate the binary two body problem instead of 
taking analytical solutions in order to minimize the influence of 
the numerical integration algorithm on the results. 
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For the non-resonant orbits, we chose the semimajor 
axis of the planet to have the values of 0.2 AU and 0.3 
AU. The eccentricity of the planet's orbit in all cases was 
chosen to be or 0.15, and its angular variables (inclination, 
mean-anomaly, longitude of the ascending node, argument 
of pericenter) were set to zero. 

Table 3 shows the results of the calculation of ETVs for 
a circular binary. Even though MMRs did not have too much 
influence on the stability of the system, the most prominent 
ETVs were found at 3:1 and 4:1 resonances. Figure [5] shows 
these ETVs for a 1 Jupiter-mass planet. 

In order to assess the detectability of our ETV signals 
with CoRoT and Kepler, we compared our results with the 
values of the Detectable T iming Amplitu d es (DT A) of these 
telescopes. As shown by ISvbilskv et al.l (|2010l ). given the 
typical photometric error, the value of the DTA for CoRoT 
is approximately 4 sec for a 12 mag star and 16 sec for 
a star with 15.5 mag. Similar analysis suggests that stars 
with magnitudes between 9 and 14.5 mag will have DTAs 
ranging from 0.5 to 4 sec with Kepler. 

Figure|6]shows a comparison between the values of DTA 
calculated by S ybilsky ct al. (2 01(]i ) and the amplitudes of 
ETVs obtained from our simulations. The x-axes in this fig- 
ure represent the planet's initial semimajor axis and the y- 
axes show the maximum amplitude of ETVs. In cases where 
the value of ETV is larger than 40 sec, this value is indicated 
on the top of its corresponding bar. As a point of compar- 
ison, the maximum values of the DTAs for CoRoT (16 sec 
for a 15.5 mag star) and Kepler (4 sec for a 14.5 mag star) 
are also shown. When the ETV signal is higher than these 
observational thresholds, we assume that a planet may be 
indirectly detectable via eclipse timing of the secondary star. 

As shown in this figure, planets in mean-motion res- 
onances present the highest probability of detection. This 
is expected as these planets create the largest variations in 
eclipse timing. Also, as indicated by the left column, binary 
model 1 with low-mass stars present the best prospect for 
detecting P-type Jovian planets. As the masses of the binary 
stars increase, the prospect of detection shifts toward larger 
planets. For instance, in binary model 3 (right column), a 5 
Jupiter-mass planet can be detected in orbits smaller than 
0.2 AU. As expected, the highest prospect of detection is 
for a 10 Jupiter-mass object (lower right panel). Figure 6 
also shows that while planets in circular orbits produce the 
largest ETVs when in resonance, slight eccentricity in the 
orbit of the planet increases the prospect of its detection 
and extends its detectability to large distances. The latter 
is more pronounced for larger planets around more massive 
binaries as in these cases, the planetary orbit stays stable 
when its orbital eccentricity is slightly increased. 



5 DISCUSSION 



Table 3. Systems with ETVs that are potentially detectable with 
CoRoT and Kepler. 



Model 


^planet 


a 


e 


ETV 




(Mj) 


(AU) 




(sec) 


model 1 


1 


0.104 





316 




1 


0.126 


0, 0.15 


40, 69 




5 


0.126 


0, 0.15 


179, 349 




5 


0.2 


0, 0.15 


25, 51 




5 


0.3 


0.15 


28 




10 


0.126 


0,0.15 


339, 730 




10 


0.2 


0, 0.15 


46, 97 




10 


0.3 


0.15 


54 


model 2 


1 


0.126 


0.15 


32 




5 


0.126 


0, 0.15 


55, 161 




10 


0.126 





108 




10 


0.2 


0, 0.15 


17, 21 


model 3 


1 


0.104 





51 




5 


0.126 


0, 0.15 


31, 59 




10 


0.126 


0, 0.15 


60, 117 




10 


0.2 


0.15 


16 




Eclipse number 




Eclipse number 

Figure 5. The graphs of ETVs for a 1 Mj planet in a circular 
orbit in the binary model 3. The planet is in a 3:1 MMR in the 
top panel and a 4:1 MMR in the bottom one. 



We have studied the prospects of the detection of circumbi- 
nary planets with CoRoT and Kepler space telescopes using 
the variations of binary eclipse timing. The uninterrupted 
high precision photometr y of these telescopes (|Koch et al.l 
120041 :1X1 onso et al. 1 12008| ) has given them unique capabil- 
ity for detecting small variations in eclipse timing measure- 
ments. We have calculated such variations in several bi- 
nary models with circumbinary planets and compared the 



results with the temporal sensitivity of these telescopes. 
As expected, the prospect of detection is higher for plan- 
ets in resonant orbits as these objects create larger ETVs. 
This is more pronounced when the stars of the binary have 
low masses. This result is consistent with the findings by 
[Schneider fc Dovld (|l995l ). However, it is important to note 
that low-mass binaries may be hard to identify. 

Our study indicates that around solar-mass binaries. 
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Figure 6. Graphs of the maximum ampHtude of ETVs for all three binary models with e^in = 0. The locations of bars on the horizontal 
axis correspond to 3:1 MMR at 0.104 AU, 4:1 MMR at 0.126 AU, and two non-resonant orbits at 0.2 and 0.3 AU. Graphs are shown 
for two values of the planet' s orbital eccentricity (0 and 0.15). As a point of comparison, the maximum values of DTAs for CoRoT and 
Kepler llSvbilskv et al.ll2010l'l are also shown. 



planets with sizes down to 1 Jupiter-mass can produce de- 
tectable signals. In this case, planets have to move in almost 
circular orbit close (a ~ 2abin) to binary's center of mass in 
order to maintain stability. More massive planets on higher 
eccentricities may orbit the binary at distances up to ap- 
proximately 6abin and still be able to produce ETVs that 
are detectable by Kepler. 



Result of our simulations indicate that although the or- 
bital stability of the planet is strongly affected by increasing 
its eccentricity and the eccentricity of the binary, slight de- 
viations from circular orbits, in particular in the orbit of 
the planet, result in its periodic close approaches to the bi- 
nary and creating large ETVs. Our study suggests that when 
not in a resonance, in general, giant planets on slightly ec- 
centric circumbinary orbits show bigger prospects for hav- 
ing detectable ETVs. With their high precision photometry 
and long duration of operation, CoRoT and Kepler are well 
suited for indirect planet detection via ETVs, and have the 
capability of detecting such planets within the durations of 
their operation. 
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